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SUMMARY 


Using a North American Navion as a flight test vehicle, a study 
was made of the practicality of experimentally determining all longitu- 
dina) stability derivatives uniquely by releasing a large weight in flight. 
The solutions for the stability derivatives were attempted by matching 
the transient response of the aircraft with solutions to the equations of 


motion generated on an analog computer. 


Due to instrumentation difficulties definite verification of the 
practicality of this procedure was not achieved. However, the experi- 
mental results indicated that the stability derivatives could be found 


uniquely by releasing a large weight in flight. 


Unique longitudinal stability derivatives were found from simu- 
lated data generated on the analog computer. Better accuracy for these 
derivatives was achieved if the transient response to an elevator step 


was analysed in addition to the weight release response. 


It is considered that this procedure for the unique determination 


of longitudinal stability derivatives is suitable for further investigation. 
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INTRODUCTION 


In recent years, dynamic flight testing has received progres- 
sively increasing emphasis. The objectives of dynamic testing are 
several in number. However, receiving the greatest attention is the 
determination of stability derivatives from the transient response of an 
aircraft to an applied forcing function. There are many advantages to 
this type of flight testing. Among the more important advantages are 
the short flight time and ble tiwaiy simple flight techniques required. 


In many cases, dynamic flight testing offers the only means of finding 


certain stability derivatives. 


The basic approach in using dynamic flight testing to determine 
stability derivatives is simply described. <A time history of the air- 
craft's transient response to some input forcing function is recorded. 

It is assumed that the form of a set of linear, constant coefficient, dif- 
ferential equations describing the system is known. The time history 


is then used to evaluate the constant coefficients in these equations. 


A basic problem that is present in dynamic flight testing is the 
unique determination of all longitudinal stability derivatives. The most 
common procedure has been to assume some theoretical relation be- 
tween these derivatives. This report is concerned with the unique de- 


termination of all longitudinal stability derivatives without any such 





assumptions. This is accomplished by using a weight release a6 a 
forcing function to produce the aircraft's transient response. This 
relieves any linear dependency between the various longitudinal stability 


derivatives. 


This investigation was carried out at the James Forrestal 
Research Center, Princeton University, Princeton, New Jersey, during 


the 1958-1959 academic year. 


*~ 
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EQUIPMENT 


Airborne Instrumentation 

The test aircraft employed was a North American Navion, 
which is a four place, all metal, low wing aircraft. The Navion is 
powered by a Continental E-185 engine rated at 185 HP at 2300 RPM 
and 29 inches manifold pressure. <A three view drawing of the aircraft 
is shown in Fig. 1, and the general specifications are recorded in 


Table I. 


The aircraft was modified by the installation of a Navy Mark 
51 bomb rack which was placed on the lower center line of the fuselage, 
aft of the nose gear wheel well. The bomb rack was installed so that 
the release mechanism could be actuated either electrically by a sole- 
noid or manually by a cable. A 566.5 pound solid steel bar was equipped 
to hang from the bomb rack. This stéel bar was 9.5 inches in diameter 
and 27 inches long. Henceforth in this report, the steel bar will be re- 
ferred to as the bomb. Figs. 2 and 3 show the bomb installed on the 


aircraft. 


Standard aircraft instruments were used to determine altitude 
and airspeed. In addition, the aircraft was instrumented to measure 
normal acceleration, pitch rate and angle of attack. These variables 
weré fed to a telemetering system and recorded at a ground station. 


Also, the instant of bomb release was recorded. Fig. 6 is a schematic 


wiring diagram for the instrumentation system. 


The accelerometer used was a Giannini Model 251)3-1-20 with 
atrange of +1g. The sensitive axis of this instrument was canted 
laterally at 45°. With this arrangement, the accelerometer will only 
respond toa 45° component of normal acceleration. This effectively 
multiplies the range of the instrument by the square root of two for 
normal accelerations. The full scale range is therefore = 1a g's. 

An initial instaNation with a z3 g's Genesco accelerometer was made, 
but unsatisfactory data was obtained. No other small range accelerome- 
ter was available which could be employed with the telemetering equip - 
ment. A 2} g accelerometer canted laterally at 45° will of course 
respond to any lateral accelerations. However, the instrument provided 
reasonably good results. Full scale voltage for the accelerometer poten- 
tiometer was provided from the 28 volt battery bus. The resistance of 
the potentiometer is 2000 ohms. A bias battery was included to adjust 


the steady state voltage to the desired level. 


The pitch rate was measured with a Honeywell gyro Air Force 
Part No. JG 7005 A-24. This type instrument operates with a rotor 
power supply of single phase, 400 cycle, 115 volt alternating current. 
The transducer is of the potentiometer type with a resistance of 530 
ohms. The axis of the gyro was carefully checked to align with the air- 
craft axis. The potentiometer was supplied with five volts from the 
battery bus. These same five volts were supplied to the telemeter equip- 


ment as a full scale reference. 





The angle of attack vane was mounted on a four foot boom which 
was located on the left wing tip. <A balsa vane with a direct connection to 
a helipotentiometer was employed to reduce the moment of inertia and 
friction of the system. The potentiometer was supplied with 28 volts and 


a bias battery. 


The telemeter unit was the pulse width type and was built by the 
Applied Science Corporation of Princeton. This unit contains 43 channels 
and a normal sampling rate of 20 per second on each channel. Three 


channels were used for angle of attack and pitch rate while six channels 


were used for the normal acceleration. 


The telemetering unit is designed to use five volts as a full scale 
reference voltage. When properly instrumented, these same five volts 
should be the full scale reference for each transducer. Small variation of 
the aircraft voltage will not effect instrument calibrations with this 
arrangement. However, it was desired to improve the sensitivity of the 
angle of attack and normal acceleration measurements; therefore, 28 
volts were employed as full scale. Since the telemeter will accept only 
zero to five volt signals, a bias battery was employed to set the steady 
state potentiometer inputs within this range. This has the disadvantage 
that a change in the 28 volt bus will change the steady state output of the 
angle of attack vane and the normal accelerometer. However, the sensi- 


tivity will not be effected. 





Ground Station Equipment 


The ground station equipment included a telemeter receiver, a 
magnetic tape recorder, an analog computer and a Sanborn model record- 


er. These components are shown in Figs. 4 and 5. 


The telemeter receiver was built by the Applied Science Corpor- 
ation of Princeton. This unit is designed to receive the coded output of 
the telemeter transmitter, decode the signals, mix channels as desired 


and transmit the data to a recording device. 


An Ampex Model 309C dual track tape recorder was employed to 
record the data. This provided a convenient means of retaining data until 
ready for further analysis. The fact that the data was in an electrical 
form was a great advantage when the use of an analog computer is contem- 
plated. All voice transmissions were also recorded on the second tape 


channel. 


A Goodyear Aircraft Corporation Model L 3 (GEDA) linear 
electronic differential analyser was employed to analyse the data. This 
analog computer provides twenty-four stabilized DC amplifiers with an 
open loop gain of 5x 10’. The computer had a guaranteed accuracy of 
one percent. Provisions for accurately setting the computer board poten- 


tiometers was available. 


A Sanborn Model 154-100 B four channel recorder was used 
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for visual presentation of the data. This instrument had a high natural 


frequency, and the low frequency response was flat to zero. 


Instrument Calibration 


Three components of the instrumentation system required cali- 
bration. These were the normal accelerometer, pitch rate gyro and angle 
of attack vane. The aircraft airspeed indication and altimeter were not 
calibrated as previous investigation indicated small errors in these 


systems. 


~ 
= 


The accelerometer was calibrated by tilting it through J 90°. 
The angle of tilt and the percent of full scale voltage were recorded. 
This checked very closely with the factory calibration which was there- 
fore accepted as valid. Fig. 7 shows the sine of the angle of tilt versus 
the potentiometer output in percent full scale voltage. The slope of this 
curve was found to be 2.115 g's for full scale voltage. Since the instru- 
ment was canted at 45°, this value was divided by the sine of 45° to 
give 2.995 g's for full scale. As installed in the aircraft 28 volts were 
applied across the accelerometer potentiometer with a five volt portion 
of this representing the telemetering full scale. As received at the tele- 
metering ground station, full scale was therefore: 5/28 x 2.995 = 


534 g's. 


Part way through the flight test program, it was discovered 
that when releasing the bomb electrically, the bomb release solenoid 
reduced the battery bus voltage by .6 volts. This reduced voltage lasted 


as long as the spring loaded bomb release switch was held down. This 





ie 


effectively reduced the accelerometer zero by 12% of telemetering full 
scale. This is indicated by the negative peak on the flight 4 accelerom- 


eter trace in Fig. 14. 


Henceforth, the bomb was released manually. However, in the 
two additional drops which were accomplished the accelerometer was 
found to have a dead zone in the steady state position. Time would not 
permit additional drops. Careful recalibration of the instrument showed 
it was linear outside the dead zone, and the zero was displaced .6 g. 
The new calibration is shown on Fig. 8 as .518 g's for telemeter full 


scale. 


The dynamic characterisitics of the accelerometer were deter- 
mined as follows. The instrument was placed on a board with the sensi- 
tive axis vertical. The support was then sharply removed, and the tran- 
sient response of the instrument recorded during its free fall. A drop of 
approximately four inches was sufficient for the instrument to reach a 
zero g steady state. The instrument dynamic characteristics were de- 
termined by matching an analog computer trace to the instrument re- 
sponse. The damping ratio of .481 and the natural frequency of 6.00 
cycles per second could then be found from the computer potentiometer 


settings. 


The pitch rate gyro was calibrated by placing it on a constant 
speed turntable and by recording the revolutions per minute and the poten- 
tiometer output. In order to decrease the full scale range, a portion at 
each end of the potentiometer was shorted with metallic paint. The final 
calibration curve is shown on Fig. 9. Full scale deflection represents 


£ 8.5° per second. 
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The gyro dynamic characteristics were determined by deflect- 
ing the rotor and by recording its return to the equilibrium condition. 
The analog computer was used to determine its damping ratio of .0597 


and its natural frequency of 4.58 cycles per second. 


The angle of attack vane was calibrated in two ways, On the 
ground, a protractor was mounted in such a way that the angular deflec- 
tion of the vane could be read. The angle and potentiometer output were 
recorded and are shown on Fig. 10. An in-flight calibration was made 
by using the telemetering receiver to record the vane position. With 
the aircraft in steady level flight, the aircraft attitude was read with a 
bubble level. By flying at varying speeds, a direct plot of angle of 
attack versus percent telemeter full scale is then available. Fig. 10 
shows both calibrations. The ground calibration indicated telemetering 
full scale was 9.04 degrees while the inflight calibration gave 7.95 de- 
grees for telemetering full scale. The dynamic characteristics of the 
instrument were not determined due to poor data obtained from thie 


instrument. 


In addition to the instrument calibration, it was necessary to 
determine the aircraft weight, c.g. location, and moment of inertia. 
Wheel scales were employed to determine the aircraft total weight. 
The aircraft was then placed on two wing jacks and supported at the tail 
with an accurate balance scale. Thec.g. location was then determined 
from the location of the support points. These measurements were 
made with the pilot in the aircraft andthe wheels raised. The test air- 


craft was equipped with a 55 pound sliding lead weight which was 
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mounted in a hollow tube in the fuselage and which could be moved 17.2 
feet through a crank arrangement in the cockpit. The flight tests were 
made with the lead weight in the aft and in a mid position. The center of 
gravity under these conditions is 29.2 and 32.5% mac. While measuring 
horizontal distances with the aircraft in the flight attitude, the bomb c.g. 
by 3.4 and 1.1 inches for the aft and mid c.g. positions respectively. 
Because the bomb c.g. is at a reasonable distance vertically from the 
aircraft c.g., these distances are dependent upon an accurate attitude 
measurement of the flight conditions at the time of the bomb drop. This 
was accomplished by a flight with the bomb installed using a bubble level. 
The aircraft was then placed in the same attitude while on jacks on the 
ground and the aircraft and bomb c.g. positions were determined. The 
moment of inertia of the aircraft was determined by oscillating the air- 
craft about the jack points with a spring mounted on the tail. Fig. 11 
shows the configuration employed. The damping under these conditions 
proved to be negligible. The natural frequency was determined by timing 
50 cycles with a stop watch. Knowing the natural frequency and the 
spring constant of the system, the moment of inertia about the jack points 
was determined. Calibration of the springs showed them to be linear for 
the load range required. The moment of inertia was then transferred to 
the aircraft c.g. The moments of inertia about the aircraft c.g. were 
found to be 3536 and 3126 slug ft.© with the c.g. in the aft and mid 


positions respectively. 
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PROCEDURE 


The procedures employed in this investigation were basically 
very simple. A steady state flight condition was established; the con- 
trols were locked; and the bomb was dropped. The initial airspeed, 
pressure altitude, and temperature were recorded by the pilot. The 
transient response of the aircraft, as sensed by the instruments, was 
transmitted through the telemetering units and recorded on magnetic 


tape. 


In detail, the procedures were carried out as follows. Since 
it was desired to recover the bomb, a water drop was impractical. In 
order to drop over land with reasonable safety, only low altitude drops 
at less than 500 feet were considered practical. In addition, a high alti- 


tude drop would require a major excavation to recover the bomb. 


Since very smooth air was known to be absolutely necessary, 
it was decided to make all drop flights very near sunrise. With gener- 
ally good weather, the air was smooth below 500 feet at approximately 
every Other day. This gave an average of about two satisfactory morn- 
ings a week to attempt the drop. Winds of three to five miles per hour 


rendered conditions unsatisfactory for the low altitude test flight. 


The initial intent was to secure good bomb drop data for three 
c.g. positions at a single airspeed. In addition, it was intended to 
record the aircraft response to a small elevator step for at least one 
set of conditions common with a bomb drop. A total of nine bomb drops 


were made. Due to instrumentation and weather problems, no com- 
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pletely satisfactory drop was accomplished. However, two drops were 
made where it was felt that the unknowns could be removed with reason- 
able confidence. Also, one elevator step was recorded that was felt to 


be usable. The data for these runs are shown in Figs. 14 through 16. 


The pre-flight checks required for a test flight were practically 
non-existent. It was found desirable to check that the aircraft and instru- 
mentation were in good condition the afternoon prior to a test flight. The 
weather predictions were checked and arrangements made to gain access 


to all necessary equipment for the drop on the following morning. 


Prior to take-off, radio contact was made with the ground sta- 
tion to ascertain proper functioning of the telemetering installation. 
After take-off, the lead weight was placed in the desired location, and 
the aircraft was trimmed for level flight. When proper power settings 
were established the ground station was advised that a run was being 
commenced. The elevator was locked approximately 30 seconds prior 
to the bomb release with a clamp placed on the control column. If air- 
speed and attitude remained constant for the final 30 seconds and no 
gusts were encountered, the bomb was released. It was found that 200 
to 400 feet was the most practical altitude range for bomb release. 
Higher altitudes made bombing accuracy uncertain. Lower altitudes 
were found to lack the necessary smooth air and required excessive 


care in maintaining safe clearances. 


Once take-off was made, the bomb was dropped regardless of 
any instrumentation or telemetering problems, This required that 
streamlining be remounted on the bomb, but it was felt that it was 


unwise to make a landing with the bomb installed. 
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DISCUSSION AND RESULTS 


Theory 

In studying Equations of Motion of the longitudinal transient 
response of an aircraft, it is conventional to assume that velocity 
remains constant. Also, the aerodynamic forces can be expressed by 
a Taylor expansion about the equilibrium condition retaining only linear 
terms. Under these conditions Ref. 1 develops the equations of motion 


as follows: 


CL 
(1) —~* a + da - dO = €; (t) 


(2) Cm @ + Cy. dat Cy. dG - hd®@O = f(t) 


a da d© 


In addition, normal acceleration can be expressed as: 


(>) te n= dO - da 

The particular forcing function in this investigation can lead 
to some confusion in considering initial conditions and in placing the 
equations in non-dimensional form. In this investigation, the aircraft 
transient response was induced by dropping a weight with the controls 
fixed. Clearly, before the weight is dropped, the lift is equal to the 
total weight of the aircraft and bomb. However, the mass of the air- 
craft during the transient response does not include the bomb weight. 
To avoid confusion, it is best to think of the aircraft configuration for 
t <0 to be that of the basic aircraft (no bomb) with its c.g. deter- 


mined, plus the effects of the bomb acting in relation to the basic air- 
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craft c.g. I€ the airplane plus bomb configuration is trimmed at t <0, 
then releasing the bomb is equivalent to applying an up force at the bomb 
c.g- Thus, one forcing function in the lift equation is an up force equal 
to the bomb weight, and in the moment equation, we have a positive 
pitching moment equal to the bomb weight times the distance between 


the basic aircraft c.g. andthe bomb c. g. 


7 eo Wy = bomb weight 


Cc 2W 2 Wp 
L. s + 
zs psve eet sv2 


i] 2W 
od psv? 


It is pointed out that the initial angle of attack is prescribed 
by it. ag defined above. However, the aerodynamic time (7) and 
the dimensionless moment of inertia (h) are defined by using the air- 
craft mass during the transient response. It is impossible therefore 


to duplicate all conditions using any forcing function other than that 


described previously. 


In addition, it is necessary to examine the expression for 
normal acceleration. This is derived from the expressions developed 


in Ref. 1: 


(7) Foe nW« mV (0 -a@} 


If this expression is changed to non-dimensional form, the 


following results: 





lg 


(8) trnWwW 3 mV (dod - da) 


Op on a=. & (40 - a 
ey ( ) 


hep. zs (d@ -da) 


(10) xn 
The forcing functions in dimensional form are the bomb 
weight in the lift equation; the bomb weight times the distance between 
the aircraft and the bomb center of gravity in the moment equation. 
In non-dimensional form, ‘the forcing functions become: 


For t<to f(t) ” f,t z 0 


For t>to 
Gx} 
(11) f,) (t) = - "> ». L, 
psv- W 2 
2Ww 
(12) f, (t) . b“b «= - Wp a a 
Sv G W Cc 


Define the quantities: 


Then: f, (t) = -**L and f,(t) «= - AC 


All time histories obtained in transient responses are in 


real time, 
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Cc e \ 
ioe es a+ tf - oe oo 
2T T 
Cc. Cc ce . “ 
(14) ale ae qa a + d& 6 -h = - mM 
T J T 7? 


These expressions are identical to those obtained using an elevator 
step input except for the definitions of the forcing functions. For an 


elevator step input: 


It is pointed out again that all conditions cannot be matched for the 
two types of forcing functions. In practical term., initial angle of 
attack and velocity cannot be identical for the two cases at the same 


density altitude. 


There is an additional point of view which can be taken about 
the initial conditions of the bomb drop. Equations (1) and (2) are 
perturbation equations. By considering fs (t) and f, (t) zero before 
the bomb drop, this establishes a, and dG, as zero. Ina linear 
perturbation problem where the initial and final solutions are both 
steady state conditions, it is arbitrary as to which is considered the 
undisturbed state. It is equally valid, in the aircraft transient re- 
sponse to a bomb drop, to consider the final steady state as the undis- 
turbed condition with @ and dO both zero. The physical motion 


of the aircraft is a steady pitch rate. 
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This can be more clearly demonstrated with the use of equa- 
tions, The variables can be expressed as the sum of a steady state 


and a transient solution. 


~ >i 
a= @,. + ay dG 3s dG... + dc, 


Under the previous arrangement when t = O: 


a, x O = ~a, dco =z O d©,. = -dG, 


If the final steady state is- defined as the undisturbed condition, the 


equations become : 


oon eee.” ° asa, duo = de, 


Under these conditions f, (t) and f, (t) are zero after the bomb drop 
and are the negatives of the previous forcing functions before the bomb 
drop. Thatis, the aircraft before the bomb drop has a force and 
moment acting until release. The initial perturbation angle of attack 
and pitch rate are unknown in this case. An analytic solution of the 
equations will show that the steady state solutions in the first case 


are the negatives of the initial condition in the second case. 


The equations for the steady state solutions can be found from 
equations (1) and (2) for the original arrangement. 


ai 


C 
(15) d@,. = a/z lth - m 


f, (t) 


C 
(6) ag, = ‘2 4 Mao 18H) 
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The initial values for d@ and a@ “inthe second case are repre- 
sented by identical expressionsas above. However, the forcing 


functions change signe. 


In the real aircraft problem, there is some disadvantage in 
considering the final state as the undisturbed condition. The aircraft 
may not reach a complete steady state in the time period where veloc- 
ity is essentially unchanged. There is also a mental difficulty in 
thinking of a steady pitch rate as an undisturbed condition in relation 


to level flight. ° 


Linear Dependence 


The determination of unique values for stability derivatives 
from flight test data is a problem that arises in any experimental pro- 
gram. The stability derivatives which can be determined uniquely is 
controlled by the linear dependence of the equations describing the 
system. This subject is considered in detail in Refs. 2 and 3. For 
the longitudinal motions considered in this report, the problem can 


be resolved without great difficulty. 


In equations (13) and (14), it is considered that all flight 
variables and the forcing functions are known and are not zero. 
Under these circumstances, there is no linear dependency between 
the two equations. All stability derivatives and the dimensionless 
moment of inertia (h) can theoretically be determined. The ability 
_todetermine C_ and C_ independently by dropping a 


da d© 
weight is the basic reason this investigation was undertaken, The 





Ze 


use of the elevator as a forcing function does not accomplish this for 
one reason. The tail force introduced into the lift equation is so smail 
that its effect on the aircraft motion is negligible. Equation (13) may 
effectively be equated to zero and a relation between the variables a, 
@ and 6 results. Therefore, the variables in equation (14) are not 
linearly independent, and it is impossible to find all the stability deriv- 


atives in this equation uniquely. A complete development of what can 


be found is contained in Ref. 2. 


In using a bomb drop as a forcing function, there is an addi- 
tional advantage. The exact nature of the forcing functions are known 
as a force and a moment. With an elevator step, this is not the case. 
It is pointed out that the moment produced by the bomb release does not 


have to be known to solve the problem if h is known. 


The facts stated above can be demonstrated by the use of the 
equations as follows. If equation (13) is used to eliminate 0 from 


equation (14), the following results: 





c C C C C C 4 
(17) ( ect Ly a) a + | Mag + Mga -h_ ta CG 
é 2 T 


v 





“ T 


AC, Ge 
de co ACY -A C 





23 


wr * Cy and 4©m are known and are not zero, all quantities in 
parenthesis can be determined with a time history of a, @ and a . 
If h is known, it is not necessary to know the value of ACm to 
determine these quantities. Inspection of the quantities in parenthesis 
shows that all derivatives are determined uniguely. However, if 

4Ci, = 0, as is the case for an elevator input, the derivatives cannot 


ACL, are both known. It is also 


be found uniquely even if h and 
pointed out that ACr, can be identically zero in the case of the bomb 


drop without destroying the uniqueness of the solution. 


ee ues y Hmedtigetion 

The previous analysis verified that a solution was possible for 
all longitudinal stability derivatives by dropping a weight. However, it 
was still necessary to verify that the transient response of the aircraft 
would be of sufficient magnitude that the instrumentation available could 
sense the flight variables with the necessary accuracy. It was deter- 
mined that about 500 pounds could be carried externally on the test air - 
craft. This represented approximately 20% of the basic aircraft weight. 
It was estimated that the bomb and aircraft centers of gravity would be 
approximately eight inches apart. Refs. 4 and 5 provided reasonable 
estimates for the longitudinal stability derivatives. With this informa- 
tion, it was possible to estimate the response of the aircraft when the 


bomb was dropped. 
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Analog of the Test Aircraft 

It was decided that the use of an analog computer was the most 
convenient method to make a preliminary analysis. Equations (13) and 
(14) were set up on the computer and the time history of the aircraft 
simulated flight parameters were placed on a Sanborn recorder. In 
addition, equation (7) was set up to give normal acceleration. See Fig. 
12 for a schematic wiring diagram of the analog computer setup. This 
arrangement allowed variation of the forcing functions and all stability 
derivatives as desired. This arrangement assumed the aircraft was 


perfectly described by the equations and that the instruments would per- 


fectly record the flight variables a, 8 and n. 


Solution for Stability Derivatives from Analog Data 

Ik was found that the stability derivatives could be found with 
reasonably good accuracy using the very simple statistical methods dis- 
cussed in Ref. 6. However, it was found that greater accuracy in much 
less time could be achieved. by curve matching two computer solutions. 
One person would set in arbitrary stability derivatives and produce a 
time history of aircraft motions. A second person would then reproduce 
the curves as closely as possible and determine the stability derivatives 
from the potentiometer settings. Excellent results were obtained in this 

Cc 


manmmner for the derivative “mao . Good results were found for M vy 


while only fair accuracy could be expected for sec yer and Cg , 
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It was also discovered that the variation in angle of attack 
would be quite small. This is of course related to the difficulties in 
determining mo and “Ly - Improved results could be obtained 
by increasing AC, that is, moving the bomb further from the air- 
craft center of gravity. Unfortunately, the location of the bomb on the 
aircraft was quite closely fixed by the landing gear wheel wells. The 
remaining solution was to shift the aircraft center of gravity. This is 
not nearly as satisfactory since, in practice, shifting weight ait in the 
aircraft increases its moment of inertia and rapidly decreases Cas : 
Both of these effects tend to reduce the size of the aircraft's response. 
The analog computer demonstrated that moving the aircraft c.g. aft 
increased the response qualitatively about half as much as moving the 


bomb c.g. forward. The character of the changes were distinct how- 


ever. 


In addition, the use of elevator input to create a sizeable 
moment input was investigated. It was found in this case that the re- 
sponse was effected most strongly by Om and Sm, - It was 
known from the previous theoretical discussion that no unique solution 
existed for ae ; Ome and man: However, the computer 
demonstrated that if a reasonable approximation was available for rate 

Cc 


derivatives, good results could be obtained for My as well as 


— 
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Using both type inputs, Erne, ‘ Gin 5 and ue could be 
reproduced toa close percentage of the initial value. omg could 
not be reproduced with a high percentage of accuracy, but in terms of 
absolute value, it was well defined. This difficulty with mae is 
somewhat inherent due to its small size compared to Oman and 
Cm,. Ae Ref. 7 states, the accuracy with which a stability derivative 
can be found is determined by the magnitude of its influence on the motion 


of the aircraft. TBA fy has a relatively small influence on the air- 


craft's motions and therefore will not be found with high accuracy. 


The variables which were to be measured included nermal 
acceleration, pitch rate and angle of attack. Any two of the three will 
provide a complete solution as is shown by equations (1), (2), and (3). 
Under the conditions set into the analog computer for the bomb drop, 
the time histories indicated the following change from the steady state 
in the variables to be measured: 

Normal acceleration: 0to .4¢g 


Pitch rate: 0 to 6° per second 


Angle of attack: -.5° to +2° 


It was felt that acceleration and pitch rate could be measured 
in this range with reasonable accuracy. However, it was felt that 
angle of attack could probably not be measured with sufficient accuracy 


to give any useful information. 
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eee re eee ey aeegetian 
The following facts were indicated by the use of the analog com- 
puter in the preliminary investigation. 

a) The longitudinal stability derivatives could be found uniquely by 
dropping a weight appraximately 20% of the aircraft weight. 

b) The weight should be located further from the aircraft center of 
gravity. <A distance of 1.00 to 1.75 feet would be satisfactory for the 
Navion. (2 = .176 to .305) 

c) The derivative map would be determined most accurately. 
omga would be found with only fair percentage accuracy but good abso- 
lute accuracy. The use of an elevator step would assist in determining 
Cmy and City . However, as mentioned previously, all conditions 
of flight cannot be duplicated for the two forcing functions. 

d) It was likely that no reliable anple of attack information would 
be obtained. 

e) The best means of obtaining a solution was probably matching 


analog data with the measured response of the aircraft. 


The general technique for finding a set of stability derivatives 
by curve matching was believed to be as follows: 
a) With the best available estimates, create time histories with 
the analog computer of all flight variables which are available from the 
test flight. The step forcing functions are easily created by using the 


reference voltage of the computer. 
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b) Vary each stability derivative separately and make an overlay 
for each derivative showing its effect on the time histories. 

c) Carry out this procedure for each type of forcing function 
considered. 

d) Use the bomb drop time history and create the best match 
possible. 

e) Shift to an elevator step. Leave the derivatives Cmag and 
oman fixed, and create the best match possible. 

f) Repeat steps (d) and (e) until the best possible match is 


achieved. 


Some small allowance could be made for the difference in flight 
conditions between the two forcing functions. If identical indicated air 
speeds are used, the initial angle of attack for the elevator step will be 
approximately two degrees less than that for the bomb drop. This con- 
dition is largely relieved by the fact that the angle of attack increases 
three to four degrees during the elevator step while it shows very small 


change during the bomb drop. 


Results Using Flight Test Data 

As previously mentioned in the report, no completely satisfac - 
tory data was obtained. Two bomb drops and one elevator step were 
selected as providing usable information. Figs. 14 through 16 show 


this data. The flight conditions for these runs are recorded in Table Il. 
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When the bomb rack installation was completed, it was found 
that the aircraft and bomb centers of gravity were considerably closer 
together than had been initially estimated. This had a twofold disadvan- 
tage. As previously discussed, it reduces the size of the aircraft 
response, making the instrumentation requirements more severe. 
Also it makes the determination of the size of the moment forcing func- 
tion more difficult. The horizontal distance between the aircraft and 
bomb centers of gravity depend quite strongly on the aircraft attitude 
when the centers of gravity are very close together horizontally in level 
flight. A one degree error in aircraft attitude created approximately a 
20% error in the size of the pitching Gearedit when the aircraft c.g. 
was completely aft. The situation was even worse as the aircraft cen- 
ter of gravity moved forward. Since no instrument recorded flight atti- 
tude, it was felt that no accurate moment forcing function could be 
determined. Measurement of aircraft attitude with the bomb installed 
showed a variation too large to be ignored. It was therefore decided to 
regard the moment input as an unknown. As brought out in the theory, 
the problem still has a unique solution; however practical application 


showed it was much more difficult. 


The streamlining on the bomb was initially installed to reduce 
the aerodynamic drag. This was thought to be desirable in order to 
estimate the pitching moment more accurately. An estimate of .5 for 


the drag coefficient of the streamlined bomb was made based on its 
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cross sectional area. This estimate was taken from Ref. 1. The 
moment resulting from the bomb drag was determined to be quite small 
and therefore unimportaut since the pitching moment was not accurately 
known anyway. The streamlining was retained throughout this investi - 
gation in an a to reproduce results. I is not considered to be 
necessary for any future investigation of this nature using low speed air- 
craft. It should be pointed out that it has been assumed that there is no 


aerodynamic effect on the aircraft due to the bomb's presence. 


The matching of the time histories for two bomb drops and one 
elevator step are shown in Figs. 17 through 19. The stability deriva- 


tives obtained through these curves are recorded in Table II. 


In the analog computer solution, the dynamic characteristics 
of the accelerometer and pitch rate gyro were included. The computer 
circuit diagrams for these instruments are shown in Fig. 13. The use 
of filters to remove noise from the test data was considered. In gen- 
eral, it was found that the stray noise could be removed by eye better 
than by filtering when curve matching was employed as a method of 


solution. 


The effect of each stability derivative on the computer time 
histories is shown in Figs. 20 through 24. It can be seen that ACr, 


any and my are not clearly defined by the use of the bomb drop. 
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However, mag and “mag appear to be reasonably well defined. 
The elevator step improves the sy and “ma » but the failure to 
reproduce flight conditions prevents high accuracy. Figs. 25 through 
29 show the effect of each stability derivative on the response to an 


elevator step. 


In spite of the instrumentation difficulties previously discussed, 
the numerical results of this investigation are considered to be reason- 
ably good. In the elevator <. shown in Fig. 19, there are no ambigui- 
ties in angle of attack or pitch rate. Table IV records the derivatives 
as established in this investigation and those arrived at in Refs. 4 and 5.» 


Cc Cc 
Maa and lio established in this inves- 


The derivatives Cmag 
tigation fall within the band of previous evaluations. Conditions were 
not reproduced sufficiently to allow a good comparison for = : 
The indications are that (my as established in this report is at 
least 10% high. “mg was determined with the aircraft c.g. located 
at 29.2% and 32.5% mac. This gives a theoretical decrease of .17 in 


Cc 
the magnitude of “a. The experimental decrease was .22. 


The calculated and experimental values of SC, for the 
bomb drops agreed within 15%. In both cases, the calculated values 


were greater. See Table WI. 
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Proposed Improvements 

The basic improvements which should be made have already 
been discussed. It is an obvious fact that the instruments employed 
should be reliable and cover the desired range without modification. If 
telemetering is employed, it is a great advantage to use the telemetering 
full scale voltage as the reference voltage for all potentiometer trans - 


ducers. 


It would be of considerable assistance if aircraft attitude and 
elevator position were also measured. While neither is required, they 
both help to determine the quality of the steady state established prior to 
the bomb release. Ih addition, aircraft attitude could be used to estab- 
lish the moment forcing function. Elevator position would verify that the 
controls remained fixed during the bomb release and could be used as a 
forcing function if an elevator input was employed. This would also 
allow the determination of Cm & rather than “°m for the elevator 
input. | 

The last but most important —e has already been dis- 
cussed in some detail; the bomb should be moved further from the air- 
craft center of gravity. The importance of this change is shown in 
Figs. 30 through 33 which present analog computer time histories of 
pitch rate and normal acceleration. The moment input here corresponds 


to a distance of 1.0 fect between the aircraft and bomb centers of gravity. 
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The effect of each stability derivative on this type of response is shown. 
In these comparisons, the same numerical values were used as those 
employed with the flight test comparisons shown in Figs. 20 through 24, 
With larger distance from the bomb to the aircraft c.g., the moment 
input could be determined accurately. These curves show that each 
stability derivative has a larger effect on the transient response when a 
larger moment input is used. The relative importance of each deriva- 
tive is not greatly effected; -however omaa appears to show the 


preatest improvement. 


There is a characteristic of these curves which is of consider- 
able assistance in determining the derivatives by curve matching. 
There is a short time where Tig and Cma have no effect on the 
transient response. This is represented on the curves by the crossing 
of the variables n and dO as “my and =i GQ are varied. 

This means physically that the aircraft is passing through Zero angle of 
attack. Fig. 34 shows that angle of attack initially decreases and then 

increases through zero to a positive steady state value. In the vicinity 

of the minimum, “may has no effect on the transient response while 
oa, and Salis have no effect where a passes through zero. In 
the flight teste, the steady state angle of attack was negative and no 


zero value occurred. 


If the forcing functions are known accurately, approximate 
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values for the stability derivatives will locate the time periods where 

@ is very small and where the minimum occurs (da = 0). These 
C C 

regions may then be used todetermine ™d9 and ™da _  independ- 


ent of ”,,. and oy. 


| Cc 
The use of an elevator step to determine Ly initially is 
advantageous. Since an angle of attack vane can be used with this input, 
Cc 
time histories of n and a give bq directly. For the elevator 


step, equations (7) and (13) show that: 


Cy, 
ng a @ ov 


Vv at 





With SA known, the remaining derivatives can be determined with 
the analog computer as follows. If a reasonably good match is not avail- 
able initially, vary all stability derivatives except — to achieve a 
fair correlation throughout. Leave “mag fixed and vary Cmag 
and Cm g to match the regions where dq _ is very small (the mini- 
mum and the steady state of a). Then vary “28 to match the 


initial portions of the time histories. 


When a good match is achieved, shift to the elevator step input. 
C C C Cc 
Leave ~“dq@ and ~™do fixed and adjust “™aq and ~'q@ to 
achieve the best match for all three variables (n, @ and dO}. Return 


C C 
to the bomb input and leave Ly and My fixed and achieve the 









So C . 
best match possible varying "da and ~“™dd. Shift between the 


two forcing functions using these techniques until the best over-all 








match is achieved. 


While the bomb drop establishes all derivatives uniquely, the 

use of the elevator step will considerably improve the accuracy with 
ome Cr ¥ 

which a and @ can be determined. The accuracy of the 
rate derivatives is improved only slightly. Since an angle of attack | 


vane can be used successfully with an elevator step, an extra cross 


check is available. 


















CONCLUSIONS 


The proposed method of uniquely determining longitudinal 
stability derivatives for slow speed aircraft is believed to be feasible. 
Due to instrumentation difficulties, positive verification of this is lack- 
ing. However, the numerical results for the Navion contained in Table 


II are considered reasonably valid. 


The stability derivatives were determined uniquely from simu- 
lated data generated on the analop computer. The simplest and most 


accurate resuits were obtained by a curve matchings technique. 


The use of an elevator step function in addition to a weight 
release will greatly assist in the determination of the longitudinal sta- 
. bility derivatives. No consideration was given to any other types of 


forcing functions. 


The accuracy with which each stability derivative is determined 


is associated with its influence on the aircraft response. The derivative 


Omar 
most accurately determined by the release of a weight is do 


Mga is uniquely found with fair accuracy by a bomb release. 


C C K e — 
ee and Le can be found with better accuracy using an elevator 


C 
step and previously determined values for ‘40 and ‘der 













RECOMMENDATIONS 


As a result of this investigation, the following recommenda-. 


tions are proposed: 











1) An investigation of this nature should be continued to obtain 


improved experimental results. The configuration and the velocity of 


the aircraft should be varied in any additional investigation. 


2) A normal accelerometer with a full scale range of approxi- 
mately one g_ should be installed in an aircraft for longitudinal 


transient flight testing. 


3) The aircraft should be instrumented to measure elevator 





position and attitude to properly analyse the response from a bomb 


release. 


4) In this type of investigation in a Navion aircraft, a 500 pound 
weight should be at least one foot horizontally from the aircraft center 


of gravity. 
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TABLE 1 






Physical Characteristics of the Navion 







A. Wing 







: . |. Total Wing Arca (includes flaps, ailerons a 
and 19.8 it@ of fuselage) 1984.2 £t© 


&. Span 33.30m%s 
3. Mean Aerodynamic Chord 


4. Angle of Incidence 


Root 
Tip 


5. Twist 
Aerodyaamic 
Geometric 39 
6. Airfoil Section 
Root NACA 44° SR 
Tip NACA 6410R 
7. Aspect Ratio 
8. Taper Ratio 
9. Dihedral 
10. Chord 
Root 
Tip 
Gb. Vertical Tail 
1. Total Area 


2. Stabilizer Area 


3. Rudder Area 





Airfoil Section 
Root 
Tip 


5. Angie of Stabilizer with fuselage 


Horizontal Tail 


l. Total Area 


4. otabiliger Area. 
3. Elevator Area 


4, Span 
5. Mean Aerodynamic Chord 


6. Angle of Incidence 
7. Airfoil Section 
8. Aspect Ratio 


9. Taper Ratio 


General 

1. Length Overall 
2. Tail Length 

3. Weight (full fuel) 


4. Engine 


TABLE I (Cont.) 










NACA 0013.2 Mod. 
NACA 0012-64 Mod. 










2° Nose left 


~~ 


43706 st” 
28.95 £2 
14.10 £t® 
13. Vee £2 
3.34 £12 
40 
NACA 0012 
4.02 


.67 


PA es) 
16.88 ft 
2430 pounds 
Continental E-185 


185 HP at 2300 RPM, 
29’ Hg M.P. at sea level 








Date 
es spe ed mph 
Temp 


Pressure 
Altitude (ft) 


Cm. 
V; (fps) 


V (fps) 





Flight 4 


23 April 59 


94 
rae om 
550 
29.2% mac 
137.8 
136.6 


1.018 


2540 
748 


611 


. 1360 
‘3132 
. 0787 


1.296 


Test Flight Conditions 


—, 


Flight 9 





2 May 59 


96 


OG. 


150 
32.5% mac 

140.5 

139.4 

|; Os 


2538 
.720 


. 589 


Be og Wa 
3541 
0887 


1,275 


@ May 59 











94 


9° C, 


450 
32.5% mac 
137.8 
137.3 
1.006 


2525 


.607 


3541 
0887 


1.298 





Stability Derivatives from Flight Tests 


Flight 4 Flight 9 and 9B 
Date 23 April 59 @ May 59 


c.g. location 29.2% mac 32.5% mac 
oi 4.98 4.98 

oma - 825° | 260 

Cmaq Bote .073 


Cmag 167 66 


om . 00244 00463 (Flight 9, bomb drop! 


Calculated “@m 0029 0065 
for Bomb Drop 
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Fig. 2 Bomb Rack Installation 





Fag. 3 Test Aircraft 








Fig. 4 Tape Recorder and Telemeter Receiver 





Fig. 5 Analog Computer and Sanborn Recorder 
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